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Abstract. An unusual symmetry of the equation for the generating function of trans- 
verse velocity differences Av = v(x + r) — v (x) is used to obtain a closed equation for the 
probability density function P(Av,r) in strong three-dimensional turbulence. It is shown 
that the terms, mixing longitudinal and transverse components of velocity field, dominate 
the pressure contributions. The dissipation terms are closed on qualitative grounds. The 
resulting equation gives the shape of the pdf, anomalous scaling exponents and amplitudes 
of the moments of transverse velocity differences. 
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Intermittency of strong turbulence seems to be a well-established experimental fact. Man- 
ifested as anomalous scaling of velocity structure functions, this phenomenon resisted a 
theoretical description for almost three decades ( for an extensive review see [l]-[3]). Early 
attempts to attack the problem (cascade models) were based on the so called refined simi- 
larity hypothesis (RSH), connecting dissipation fluctuations with those of the differences of 
longitudinal components of velocity field. Only recently the properties of tranverse compo- 
nents of velocity field gained some interest [4] -[5]. In this paper we explore a remarkable 
feature of the equations for the generating function of tranverse components of velocity field 
discovered in [6] to derive the equation for the probability density (pdf) of transverse velocity 
differences. The solution to this equation describes the experimentally observed shape of the 
pdf, anomalous scaling exponents and the amplitudes of the structure functions. 
We are interested in statistical properties of an incompressible fluid strirred by a random 
force in the right side of the Navier-Stokes equations, defined by the pair correlation function: 

< /,(k)/,(k') >oc - -JL2) \ d ^ f> ^ + VMt - t>) (1) 

so that (f(x + r) — f(r)) 2 oc P(l — cos(kfr)). The integral scale of the problem L « 1/kf 
and d stands for space dimensionality. Consider two points x and x' and define r = x — x'. 
Assuming that the x-axis is parallel to the displacement vector r, one can find for the 
separation r in the inertial range ([1], [7], [8]) that S3 =< (A-u) 3 >oc Pr and S\ = (Av ) 3 = 
(v(x') — v(x)) 3 = where u and v are components of velocity field parallel and perpendicular 
to the x-axis (vector r), respectively. These relations , resulting from the Navier-Stokes 
equations, are dynamic constraints on any theory of turbulence. The pumping power P = 
0(1) is a constant. In what follows we will often set P — 1, L — 1 and restore the correct 
dimensionality at the end of calculations. We also have [1], [7], [8]: r 2 ~ d d T {j- d ~ 1 S2) = (d — 
l)Sl = {d- \){Kvf and for d = 3: 

6S 3t = 6Au(Av) 2 = d r rS 3 oc Pr (2) 
which will be important in what follows. 
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Following [9], [6] and [10] we consider the iV-point generating function: 

Z =< e x '- v ^ > (3) 

where the vectors x ; define the positions of the points denoted by 1 < i < N and summation 
over Xj Xj is assumed. Using the Navier-Stokes equations and incompressibility condition, 
the equation for d t Z = Aj M < (d t v^(xi))exp(\ ■ v) > can be written [6], [10]: 

dZ & Z r r 

m + ox-8x-r If + Ip + D (4) 

with the forcing, pressure and dissipation contributions //; I p ; D defined below. In what 
follows we will be mainly interested in the probability density function of the two-point 
velocity differences which is obtained from (3)-(4), setting Ai + A 2 = (see Refs. [6], [9]- 
[10]), so that Z =< exp(\ ■ U) >, where U = u(x') — u(x) = Au. Assuming , as in [9]-[10], 
that if \U\ << u rms and r << L, the inertial range variable Am is independent of the 
"large-scale" variable U+ = u(x) + u(x + r) gives: 

m + ox^r If + Ip + D (5) 

In a statistically isotropic and homogeneous flow the generating function can depend only 
on three variables: 



A • r 

Vi=r; i] 2 = = \cos(9); r] 3 = ^ \ 2 - r]l, 

In these variables {rji = r > 0): 

Zt + [d m d m + d -^d m + ^d m d m + il ^^d m - ^dl\Z = I f + I p + D (6) 

where 

I P = \< (d 2 ,iP( 2 ) ~ «9 M P(l))e A ' U > (7) 
I f = {r,l + r l l)P{l-cos{k f r))Z (8) 

D = v\» < (4,*v(2) - dlv tlI (l))e X - U >= -(vl + Vl) < (5(2) + ^(l))^ A ^ Av > (9) 
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where, to simplify notation we set <9j iQ = and v(i) = v (xi) and A x = ^2 and A y = 773. 
The dissipation rate £{x) = v(diVj) 2 . The generating function, depending on a single angle 
9, can be written as: 



so that any correlation function 

< (A«)»(At,)™ >= (d - ^ — —Z^ = V3 = 0) (11) 

Differentiating the equation (6) and setting both ^2 = ^3 = gives all known kinematic and 
dynamic constraints, outlined above. It is clear from the symmetries of the problem that the 
transverse velocity difference probability density is symmetric, i.e. P(Av,r) = P(—Av,r). 
We are interested in the equation (6)- (9) in the limit i]2 ^ 0. Let us first discuss some of 
the general properties of incompressible turbulence. Consider the "one-component forcing 
fucntion" f(x, y, z) = (f x (x, y, z), 0, 0) satisfying the incompressibility constraint (V • f = 0). 
In this case the equation (6 ) is: 

[d m d m + 2 -d V2 + - - V j(l - cos(k f r))]Z = I P + D (12) 

Then, setting r/2 = removes all information about the forcing function from the equation 
of motion. The equation (12) explicitely assumes that the flow at the scale r << L is 
statistically isotropic and homogeneous This can happen due to the pressure terms Ap = 
—ViVjViVj, effectively mixing various components of the velocity field. This universality, 
i.e. independence of the small-scale turbulence on the symmetries of the large-scale forcing 
is an assumption of this work. 

One property of equation (12) deserves discussion [6]. Neglecting for a time being I p and D 
we have in the limit ^2^0 

[d m d V2 + ~ r d m + Jq^}Z(V2 = 0,Vs,r) = (13) 

and 
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ft + - + ^^Z(0,ife,r) = tf(ife) (14) 
r r c/773 

where an arbitrary function ^(i] 3 ) can be chosen to satisfy various dynamic constraints. 
The fact that equation (14) contains only one derivative means that the underlying dynamic 
equation is linear, provided both I p and D contributions involve d m and the total order of the 
original equation n < 2. It will become clear below that this is the case in two-dimensional 
turbulence, while the situation in 3D is more complex due to dissipation contributions, absent 
in 2D. 

First, we will evaluate the pressure contribution I p which when rj 2 ^ 0, can be rewritten as: 
I P ~ Vs < (9 y p(0) - dy'p(r))exp(y/d — lr] 3 Av + r] 2 Au) > (15) 

d y p(0)-d y/ p(r) = J A;,(l-e^ r )[S M ( g ) M (A;-g) + ^( g ) u ( A; _ g ) + ^ M ( g ) w ( A; _ g )] rf 2^ g 

(16) 

and the exponent is expressed simply as: 

e V3=ivsAv+v2A V = exp [^dZTi V3 J {i- e lQxr )v{Q)d 2 Q + r] 2 J (1 - e iQ * r )u{Q)d 2 Q] (17) 

The ^-integration is carried out over the interval (— oo, oo), so that only even powers of k y 
can contribute to the integral. The most interesting feature of expressions (15)-(17) is that 
the additional /c^-factors can appear only from various correlation functions resulting from 
the expansion of the exponents. This is a consequence of the fact that here we are interested 
in the ^/-components of the pressure gradients which are perpendicular to the displacement 
vector r. 

Let us consider the N = n + m-rank tensor 

T N =< vn(x)v j2 (x) ■ ■ ■ v an (x)v ayn+ i(x + r) • • • v^^+^x + r) > (18) 

In a statistically isotropic flow this tensor can be expressed only in terms of Kronekker 
symbols Sij, componnets of the displacement vector r a and some functions A(|r|). Thus the 
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tensor can be represented as a sum of terms having a general structure A(r)5ij ■ ■ ■ 5 pc rj ■ ■ ■ r a . 
The product of h Kronekkers contains an even (2h) number of indices, while the product of 
N — 2h components of the displacement vector are the ones we are interested in. We have 

T N = J < v i (q 1 )v j (q 2 ) ■ ■ ■ v a (q n )v a (q n+1 ) ■ ■ ■ v u (q n+m ) > e tQ r (19) 

where Q = q n +i + q n +2 • • • q n +m- The tensor T/v can be represented as a sum of contribu- 
tions having the following structure: t 2 h I A(Q)®Qj ' ' ' ^ a e lQ ' r , where t 2 h is a product of the 
h Kronekker symbols. As a consequence, a typical contribution to the correlation function 
< (vi(qi)vj{q 2 ) ■ ■ -v a (q n )v a (q n+1 ) ■ ■•v U) (q n+m ) > can be written as: t 2h B(\Q\)Qj ■■■Q a where 
B(Q) oc A(Q)/Q m . Now we consider each term in the right side of (16). In the limit ^2^0 
the expansion of the exponents gives only various powers of v (?)'s. The first 0(u 2 ) term 
generates two kinds of correlation functions 

< uuv 2n >oc k 2 x k 2 y n (20) 
and < uuv 2n+1 >= 0. The second term in (16) generates a typical contribution 

< vvv m >oc k™ (21) 

if m is even- and is equal to zero if m is an odd number. Thus, both first and second terms 
in the right side of (16) are equal to zero. We are left with the third contribution, mixing 
the u and f-componnets of the velocity field. This means that, neglecting possible infra-red 
divergences the estimate for the remaining term gives 

L = b^< AuAve^+^v >= b^——Z(ri 2 = 0, 773, r) (22) 
r r or) 2 or)3 

This result was first obtained for a close- to- gaussian case of two-dimensional flow [6] . Using 
some additional assumptions about the structure of the ^(f] 3 ) the gaussian expression for 
the PDF of transverse velocity differences in two-dimensional turbulence was obtained. Here 
we have proved that, in general, the only contribution to the I p - term comes from the one 
mixing u and v components of velocity field. 
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The difficulty of the three-dimensional case comes from the dissipation contributions. This 
is easily understood on the basis of the Navier-Stokes equations. Introduce the local value 
of the kinetic energy K = v 2 /2 and Q = K(x + r) + K(r). It is clear that 



£(x+r)+£(x) ~ d t Q+d2,iV i (x+r)K(x+r)+d lii v i (x)K(x)+d 2 ,iV i (x+r)p(x+r)+d lii v i (x)p(x) 

(23) 

where the 0(v) terms are not written down. One can see that the 0(Vt> 3 ) dissipation 
contributions are much more complex than the 0(Vv 2 ) pressure terms, derived above. Thus, 

< (£(2)+£(l))e A ' Av 

^- < (v l (2)K(2)-v l (l)K(l))e x ^ v > -A, < (v i (2)K(2)s ij (2)-v i (l)K(l)s ij (l))e x ^ > +R 

where Sij(m) = and R denotes the remaining contributions coming from (23) and 

the viscous terms. Assuming that the integral scale L can enter the equations through the 
velocity differences only we make a proposition: in the limit i]2 ^ 0: 



< (8(2) + £(l))e AAv >^< AuAv^e x -* v >^< AuAv^e^ Au+v ^^ Av > (24) 

or or 

Assumption (24) is essential for the present theory. In its spirit it resembles Kolmogorov's 
refined similarity hypothesis, connecting the dissipation rate, averaged over a region of radius 
r, with (Am) 3 . It will become clear below that in the limit r\2 — > the dimensionally correct 
expression 

D° w < Av i (Av j ) 2 e r ' 2Au+V ^ Ir ' 3Av >oc r] 2 2 d r d m d 2 m Z 

or \ 

cannot appear in the equation for Z, since it leads to exponents of the structure functions 
violating the Holder inequality. It means that the model (24) can be valid only as a result 
of mutual compensation of various contributions to the expression for D defined by (9) and 
(23). At the present time we cannot prove that it is so. With the assumption (24) we have: 



7 



D oc r) 3 d m d m d r Z (25) 

Thus, in the limit r] 2 — > 0: 



[dA + 2 -d m + ^^A- - ^ A, + 6^^- + ^^^Z^ = 0, % , r) = (26) 
where unknown coefficients b and c will be determined below. Integrating (26) over r] 2 gives 

l d m + - + + + cn 3 d m d r ]Z(v2 = 0, V3 ,r) = tffo.r) (27) 

We must choose ^(773) in such a way that the generating function Z(0, 0,r) = 1. Inverse 
Laplace transform gives the resulting equation for the pdf P(Av,r) 

Since S3 = 0, the coefficients in the equation for the pdf P(V,r) = P(—V,r) (28) are chosen 



to give s 3 = \Av 1 3 = a 3 Pr with an undetermined amplitude a 3 . This is an assumption of 
the present theory, not based on rigorous theoretical considerations. Seeking the solution in 
a form S l n =< (Av) n >oc r^ n we obtain 

1 + 3/5 1.15 
£ n = _ Z—n w _ - n 29 

^ 3(l + /3n) 3(l + 0.05n) V ; 

which was derived in [10] together with (5 ~ 0.05. The equation (29) gives: P(0,r) oc r~ K 
where k = 3 ^_ 3 ^ ~ 0.4 for /3 = 0.05. Very often the experimental data are presented as 
P(X, r) where X = V/r" with 2/i = £ 2 ~ 0.696 for f3 = 0.05. This gives 

P(X = 0, r) oc r- K+ ^ « r" ' 052 (30) 

The experimental data by Sreenivassan [11] give —k + /i w —0.06. Introducing P(V, r) = 
r~ K F(V, r) leads to: 

<9F OF d 2 F 

^-^ + KV sv^ Vr ava; = (31 > 
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Next, we define Y = V/r K so that: 

BF BF d 2 F d 2 F 

(! " flr W + W*Y W + (3kY 2 — - PkYt— = (32) 

Denning — oo < y = Ln(Y) < oo, substituting this variable into (32) and evaluating the 
Fourier transform of the resulting equation gives: 

BF BF 
(l-(3)r— + (3K(ik-k 2 )F-ik(3r— = (33) 
Br Br 

with the result: F oc r 1 ^, where 

•rf*> = /»«Tq£&W> (34) 

with r/L << 1. We have to evaluate the inverse Fourier transform: 



dke~ lky e^ k) (35) 

in the limit y = 0(1) and r — > so that Ln(r/L) — > — oo. The integral can be evaluated 
eaxctly. However, the resulting expression is very involved. Expanding the denominator in 
(35) gives : 

F= / rffce 1-/3 } e ^ (36) 

J — oo 

and 

jr 1 exp(-« (37) 
with £ = V/r^e and fi(r) = A(3Kj^\Ln{r / L)\. 

To understand the range of validity of this expression, let us evaluate < V n > using the 
expression (37) for the pdf. Simple integration, neglecting 0(/3 2 ) contributions, gives: < 
yn >QC r a n w ^ ^ — (i _)_ 3/5) (77, - /3(n 2 + 2))/3. Comparing this relation with the exact 
result (30) we conclude that the expression for the pdf, calculated above, is valid in the range 
n » 1 and f3n « 1. The properties of the pdf in the range 3 < £ < 15 are demonstrated 
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on Figs. 1 and Fig. 2 for r/L = 0.1; 0.01; 0.001. Fig. 2 shows the dependence of the pdf on r 
for two values of £ = 10; 15. In this range of parameter variation the curves are reasonably 
well approximated by the the function 

P(£,r) cxe" ^ 9 (38) 

wwith a(r) r e wwith 9 m 0.12 — 0.15. All quantitative predictions (29), (30) and (38) can 
be easily veryfied exprerimentally. Indirect conformation of this result can be found in the 
experimental data by Gagne [12] on the asymmetric pdf of longitudinal velocity differences: 
the curves were approximated by an exponential with width a(r)(r) oc r e ° with 6° 0.15. 
The log-normal distribution (37), derived from (26), (28), is valid in a certain (wide but lim- 
ited) range of V- variation. It is clear from (26)-(28) that neglecting the dissipation terms 
(c oc (3 = 0) leads to = n/3, i. e. disappearence of anomalous scaling of moments of 
velocity differences. This result agrees with the well-developed phenomenology, attributing 
intermittency to the dissipation rate fluctuations: the stronger the fluctuations, the smaller 
fraction of the total space they occupy [1], [8]. To the best of our knowledge, this is the 
first work leading to multifractal distribution of velocity differences as a result of approx- 
imations made directly on the Navier-Stokes equations. The expression (37) is similar to 
the one obtained in the groundbreaking paper by Polyakov on the scale-invariance of strong 
interactions, where the multifractal scaling and the pdf were analytically derived for the first 
time [13]- [14]. In the review paper [14] Polyakov noticed that the exact result can be simply 
reproduced considering a cascade process with a heavy stream (particle) transformed into 
lighter streams at each step of the cascade (fission). Due to relativistic effects the higher the 
energy of the particle, the smaller the angle of a cone, accessible to the fragments formed as 
a result of fission. Thus, the larger the number of a cascade step, the smaller is the fraction 
of space occupied by the particles [14] . It is remarkable that the qualitatively similar picture 
was so successfully applied by Parisi and Frisch to the undertanding of scale invariance of 
strong turbulence [15]. 

To conclude: theory of turbulence needs a closure for the pressure gradient and dissipation 
terms in the equations for the probability density. The task is greatly simplified by an 
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unusual symmetry of the equations for the pdf of transverse velocity differences, resulting 
in reduction of the order of the equation for the generating function Z. It means that, 
not being able to evaluate the entire probability density P(Au, Av, r), we can derive an 
expression for P(Av,r) = /f^ P(Au, Av , r)dAu. It has been shown in this work that the 
dominant contributions to the expression for the pressure gradients comes from the terms, 
"mixing" transverse and longitudinal components of the velocity field. The self-interaction 
does not affect the dynamics of transverse velocity differences. We believe, the expression 
(22) stands on a relatively firm ground. The closure of the dissipation terms (24), (25) is 
an assumption which must be tested experimentally. The proposed closure resembles the 
kolmogorov refined similarity hypothesis connecting the dissipation rate averaged over a 
"ball" of a radius r around point x with the third power of velocity difference (Aw) 3 . The 
relation (24)- (25), though, is based on an equivalent mixed third-order structure function 
(2), untill now overlooked by studies of strong turbulence. 
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Fig.l. Ln(P(£,r)). From bottom to top: r/L=0.1; 0.01; 0.001, respectively. 



Ln(-Ln(P^,r))) 




Ln{r) 

Fig. 2. Ln(—Ln(P(^,r))) vs. Ln(r) for £ = 10 (top) and £ = 15. 
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